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Abstract 

We investigate the effect of general flavour mixing among squarks on the rare decays 
B — > Xg'j, B s — ► u + u~ and B s — B s mixing beyond the leading order in perturbation 
theory. We include all large tan /^-enhanced corrections whilst also taking into account 
the effects of general flavour mixing on the uncorrected quark mass matrix and SU(2) l x 
U(l)y breaking. For B s — > u + [i~ and B s — B s mixing we find that, in analogy to 
B — > Xg'j, there appears a focusing effect which can reduce the contribution due to 
the 5rr (and the 5ll) insertion by up to a factor of two at large tan/3 and fj, > 0. A 
dependence on 5lr and <5rx> that otherwise cancels to first order in the mass insertion 
approximation, is also reintroduced. Taking into account the current experimental 
bounds on AMb s and BR(i? s — > u + n~), we find that the insertions 5rl and 5rr can 
be significantly constrained compared to bounds obtained from B — > X s ^f only. 



1 Introduction 



Flavour changing neutral current (FCNC) processes provide a promising place to look for 
possible signals of physics beyond the Standard Model (SM). This is because the GIM mech- 
anism ensures that the SM contributions and additional effects due to "new physics" both 
enter at the one-loop level. As such the increasingly accurate experimental data obtained 
from dedicated 5-factories as well as collider experiments can be used to place constraints 
on masses and other parameters of a given new physics model. 

The process that provides some of the strictest constraints on physics beyond the SM is 
B — > Xg'-f. The current world average for the branching ratio is given by [1] 

BR(B -> X s7 ) cxp = (3.52 ± 0.30) x 10~ 4 E y > ^m b . 

The SM prediction for the branching ratio for the decay is based on a next-to-leading order 
(NLO) calculation that was completed in Refs. [2, 3] 1 

BR{B -> X s7 ) SM = (3.70 ± 0.30) x 10" 4 E 1 > —m b . 
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Using the recent results for the decay B — > X s l + l~ , the sign of the B — > X s j amplitude 
can also be determined [5] to be that of the SM contribution. These results allow further 
constraints to be placed on any new physics beyond the SM that feature large contributions 
from additional sources of flavour violation. 

Studies of physics beyond the SM such as supersymmetry (SUSY) have, until recently, 
typically focused on the inclusion of leading order (LO) effects (for example, see Ref. [6]). 
However, due to the increasing accuracy of experimental data and its relatively good agree- 
ment with the SM prediction it is becoming necessary to include at least the dominant effects 
that occur beyond the LO. 

Such effects have been studied, for example, in the two Higgs doublet model (2HDM) 
and the Minimal Supersymmetric Standard Model (MSSM). The 2HDM calculation was 
completed in Refs. [7, 8]. Studies of the MSSM contributions have tended to focus on various 
approximations and specific schemes. For instance, the results presented in Ref. [9] assume 
minimal flavour violation (MFV) and a particle spectrum in which the charginos and one of 
the top squarks are lighter than the rest of the superpartners. In Refs. [10, 11] the effect 
of large tan /3-enhanced beyond leading order (BLO) corrections to the 6-quark mass and 
charged Higgs coupling on the process B — > X s ^ was explored and shown to be sizable, and a 
resummation of terms proportional to a s tan (3 was performed to keep perturbation expansion 
under control. These methods were extended to include neutral Higgs contributions and 
tan /3-enhanced corrections to the tree-level CKM matrix in Ref. [12], general electroweak 

1 A somewhat more conservative estimate is given in Ref. [4]. 
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contributions and SU(2)t x U(l)y breaking effects in Ref. [13] and the effects of general 
flavour mixing (GFM) in the soft sector in Ref. [14]. 

The calculation presented in Ref. [14] was based on constructing an effective field theory 
where the supersymmetric particles are integrated out at a scale [Isusy (in a similar manner 
to Ref. [10]). The down quark tree-level (or, in the language of Ref. [14], "bare") mass matrix 
and the effective couplings were then calculated in the context of GFM. It was found that 
taking into account these effects can significantly reduce the bounds on the flavour violating 
parameters compared to purely LO calculations as a result of a "focusing effect" [14]. 

In this Letter we extend and generalise the methods presented in [14] to B s — B s mixing 
and the decay B s — > fi + fi~. Whilst the mass difference AM# s and the branching ratio 
BK(B S — * /i + /i~) have so far remained undetermined, both processes provide possible places 
to look for signals of physics beyond the SM. In particular, large regions of the MSSM in the 
regime of large tan (5 can be explored. Since £>-factories do not run at the energy required 
to produce large quantities of the B s mesons, the best experimental constraints on both 
processes are provided by collider experiments. The current experimental bound for the 
process B s — > is provided by the D0 group at the Tevatron [15] 

BR(B S -> fi + fi~) cxp < 5.0 x 10" 7 95% C.L, 

whilst the experimental bound on AMb s is [1] 

AM^ P > 14.5PS- 1 95% C.L. 

The SM contributions to these processes are known to NLO [17, 18]. However, the largest 
source of error for both quantities is due to the hadronic matrix element fs a which has to 
be determined using either lattice calculations or QCD sum rules. The values given in the 
literature for the branching ratio for the process B s — > fi + fi~ therefore tend to vary but are 
typically of the order [19] 

BR(B S -> ^ + /i")sM = (3.2 ± 1.5) x lO" 9 . (1) 

The SM prediction for AM Bs is [20] 

AM|f = 18.0i3.7ps" 1 , (2) 

where the large hadronic uncertainty has been avoided by using the experimentally measured 
value of AM| M . It has also been pointed out that in models with minimal flavour violation 
the large uncertainty associated with the branching ratio for the decays B q — > can also 

be avoided by relating it to the experimentally measured values of AM^ [21]. 

The contributions due to effects beyond the SM on the decay B s — > arise due to 

contributions to Z and Higgs penguins as well as box diagrams mediated by charginos and 
(in the GFM framework) neutralinos. The contributions due to neutral Higgs penguins in 
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particular have been the subject of intense theoretical investigation [22, 23, 24, 25, 27, 28, 12, 
13] recently due to the tan 6 /? dependence of BR(B S — > fi + fi~). At large tan/3 it is therefore 
quite possible for BR(£> S — > n + n ) to be enhanced by a few orders of magnitude compared 
to the SM value, whilst still satisfying current experimental bounds and the restrictions 
imposed by B — > X s 7. Similar contributions arise in the B s — B s system, where the double 
Higgs penguin diagram, although strictly an NLO effect, can become comparable to the LO 
contributions in the large tan/3 limit [30, 27, 13]. 

In this Letter, in addition to the effects discussed in [14], we include the contributions of 
charginos and neutralinos when calculating corrections to the bare quark mass matrix and 
corrected vertices. We also include the effects of the modified neutral Higgs vertex when 
evaluating the contributions to B — > X s '~f. For all three processes: B s — B s mixing and the 
decays B — > X s ^ and B s — > fi + fi~, we take into account all tan /3-enhanced corrections, the 
additional electroweak and SU(2) L x U(l)y breaking effects discussed in [13] and the effects 
of GFM mixing parameters on the bare mass matrix [14]. 

2 GFM and tan (3— Enhanced Corrections 

The effects of tan /3-enhanced SUSY corrections to the down quark Yukawa coupling [32] 
and the charged Higgs coupling [10, 11] have been found to be large and their inclusion 
can produce sizable deviations from purely LO calculations. As has been pointed out in 
Refs. [22, 24, 12, 13], these corrections can also affect the structure of the CKM matrix, K, 
due to the additional unitary transformations required to transform the quark fields into a 
super-CKM basis. 

To illustrate this consider the effect of integrating out the SUSY particles on the down 
quark mass matrix. The physical down quark masses (denoted m^) are given in terms of the 
uncorrected quark masses (denoted rrffi) by the relation 2 

m d = m d 0) + 5m d . (3) 

Note that, relative to Ref. [14], we have dropped the factor a s /Aix in front of 5m d because 
here we will also include chargino and neutralino corrections, in addition to the (dominant) 
SUSY QCD ones considered in Ref. [14]. 

In the physical super-CKM basis the mass matrix m d is (by definition) diagonal, but in 
general (and bmi) is not and provides a source of flavour violation. Alternatively [22, 
24, 12, 13], one can start with the "bare" super-CKM basis where mf^ is diagonal and, after 
computing the corrections, diagonalise the corrected mass matrix, which amounts to rotating 
to the physical super-CKM basis. In this approach flavour violation is introduced via unitary 

2 We will generally follow the language and conventions of Ref. [14] . 
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rotation matrices as they affect the CKM matrix, as well as the various other vertices present 
in the resulting effective theory. 

In the limit of MFV both approaches can be shown to be equivalent. For example, 
the 1 + e s tan (3 dependence (where e s will be defined shortly) that arises when the CKM 
matrix is modified in the approach described in Refs. [22, 24, 12, 13] is reproduced once 
the gluino contributions to a given process are taken into account. When performing MFV 
calculations the method presented in Ref. [13] is more convenient since the gluino corrections 
to a given vertex are solely to the flavour diagonal terms. However, in GFM scenarios the 
flavour violating gluino (and neutralino) contributions are evaluated anyway. Additionally 
the iterative procedure described in Ref. [14] is more suited to calculations where the squark 
mass matrix is diagonalised numerically. 

Here we follow the procedure described in Ref. [14]. The bare mass matrix and 
the corrections to the electroweak vertices are calculated in the physical super-CKM basis 
using an iterative procedure. The supersymmetric contributions to the process in question 
are then evaluated, taking into account the effects of the modified bare mass matrix, and 
evolved from Hsusy to the electroweak scale nw using the relevant six flavour anomalous 
dimension matrix. The electroweak contributions are then evaluated, using the uncorrected 
vertices when evaluating the NLO corrections and the corrected vertices when evaluating the 
LO contributions. Finally the combined supersymmetric and electroweak Wilson coefficients 
are evolved from ji w to fit using the five flavour anomalous dimension matrix and used to 
calculate the relevant observable for the process in question. 

Before presenting our numerical results it will be useful to consider the effects of including 
GFM contributions when calculating and the resulting effects on the Wilson coefficients 
relevant to B — > X s ^, B s — > fi + fi~ or B s — B s mixing. To this end we work in the mass inser- 
tion approximation (MIA) where the off-diagonal entries of the 6x6 squark mass matrix are 
treated as perturbations and flavour violation is communicated through mixed propagators 
proportional to the appropriate off-diagonal element. In our numerical analysis MIA is not 
assumed and all the squark mass matrices are diagonalised numerically. 

Departures from MFV can be measured in terms of the parameters 5ll, 8lr, &rl and 8lr 
definitions of which can be found in Ref. [14]. If one delta is varied at a time the diagonal 
terms of the bare mass matrix are given by the well known result [32] 



where i — 1,2,3 (or d, s, b), Y t is the Yukawa coupling of the top quark and the presence 
of the Kronecker 5-function 5 i3 reflects the fact that the chargino contribution proportional 
to the top quark Yukawa coupling only effects the bottom quark mass. (Corrections to the 
strange and down quark masses are suppressed by K ts K^ s and K t dK^ d , respectively, and are 




(4) 
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set equal to zero.) Finally, 

6s = "l^ 3 )^ 2 ^'^' €Y = ~i&kji H2{ - yiiR ' yiiL) ' (5) 

where a s is the strong coupling constant, C 2 (3) = 4/3 is the quadratic Casimir operator 
for SU(3), fi is the Higgs/higgsino mass parameter, rrig is the gluino mass and A t is the 
3, 3 element of the trilinear up-type soft term. The loop function H 2 can be found in the 
Appendix. Its arguments and some other quantities to appear below are defined as 



'"7l.LL.HH _ \l m d,LL\J™d,RR "'j.U.JW 

X d L , R ~ m 2 > X d RL ~ m 2 > ^2 > W 

where fn% LL , fh 2 dRR denote common values of the diagonal entries of 3 x 3 squark soft SUSY 
breaking terms for which we follow the conventions given in Ref. [14]. Whilst the diagonal 
elements of the soft terms have been assumed to be universal, it is relatively easy to include 
the effects of flavour dependence at the cost of clarity in the final expressions. 

It has been pointed out in Ref. [29] that if Sll and S RR are both non-zero large corrections 
to the bare strange and down quark masses can occur at third order in the MIA. In our 
numerical analysis we diagonalise the squark mass matrix numerically and therefore these 
corrections are automatically included. 

Taking into account flavour violating effects in the LR sector and ignoring the effects 
induced by other sources of flavour violation (including the CKM matrix) the off-diagonal 
elements of are given by (a more complete formula will be given in Ref. [33]) 

where i, j = 1, 2, 3 and 

tRL = -^C 2 (3)H 2 (x iR ,x iL ). (8) 
The effect of including GFM corrections to and the electroweak vertices can be rather 
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large. For example, in the case of B — > X s 7 the presence of non-zero (jn$^ can lead to 
large corrections to 5 X C-j$ that are otherwise suppressed by a factor of K^mi, [14]. Similarly 
( m d°' ) ) induces analogous corrections to the chargino contributions in the primed sector 
that are usually suppressed by m s Kl s . 

In the case of B s — > fi + fi~ , the gluino contributions to the Wilson coefficients of the scalar 
and pseudoscalar operators become (in the large tan (5 limit) 

$9C - ± 4:as tan 3 /3 
S ' P 3a m\m b mg K tb K^ s 



(9) 
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where m M denotes the mass of the /i-lepton, whilst the primed coefficients are given by 



6°C'a p = 



4a s jim^ tan 3 f3 



s ' p 3a m 2 A m s m~ g K tb K$ s 



+ ( m T) 33 ( S n R ) 23 x d R H ^ x d L , Xd R > x d 



(10) 



(for details of the operator basis we use see Ref. [26]). Substituting Eq. (4) and Eq. (7) into 
the above expressions yields the Wilson coefficients: 



5°C s ,p = ± 



4a s /im At tan 3 /? J e RL 



X d RL 



3a m 2 A K th Kl s | [1 + e s tan j3] m b 
1 



x d L 



[1 + e s tan (3} [1 + (e s + e Y Y t 2 ) tan (3} m~ 6 



(5 LR ) 23 H 2 (xj R ,Xj L ) 
(5 LL ) 23 H 3 (x iR , xg L ,x iL )\, (11) 



- _4a s ^ tan 3 /? j [e M + ] x~ dRh 

6 ° s ' p ~Ja~m\ K tb K* s \ [1 + (e s + e Y Y 2 ) tan /3] m s [dRL j * 3 ' ^ j 

+ [1 + (e, + eyyfl tan/?]* ^ ( ^ 23 ^(^, x, a , |, (12) 

where tha denotes the mass of the pseudoscalar Higgs boson and a denotes the electromag- 
netic coupling constant. 

Comparing the above expressions with the analysis of Ref. [27], the term in Eq. (4.15) 
proportional to the insertion 3 5ll needs to be corrected by a factor 

l + (e s + e Y Y t 2 )tanp 
1 + e s tan /? 

which reflects the additional contribution (which was omitted in Ref. [27]) obtained when 
including the effects of the insertion on the bare CKM matrix. 

Values of up to 0(1) for both 5ll and Srr are viable in some regions of parameter 
space [34] and can lead to large contributions to and Cp. Large values of 8rr, in 
particular, are motivated by SO(10) or 577(5) based solutions [35, 36] to the neutrino mass 
problem. As well as these corrections, the insertions 5lr and 5rl reappear once BLO effects 
are taken into account. At LO the insertions vanish due to an accidental cancellation between 



3 Since we are interested here in the 23 mixing elements, henceforth we shall adopt the conventional 
notation of (Slr) 2 3 a s Slr, (Srl) 2 3 a s Srl, etc. 
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the self energy and vertex corrections to the effective Higgs vertex. Including BLO effects 
however, allows the insertions to appear through their effects on the bare mass matrix Eq. (7). 
This effect is independent of the LO cancellation and originates from the additional unitary 
transformations that are required to transform between the bare super-CKM basis and the 
physical super-CKM basis. Additionally, since the corrections proportional to Slr and 8rl 
depend on rrig, rather than the strange or bottom quark mass, the enhancement by m g /mb 
can compensate for the a^-dependence of the Wilson coefficient. 

As discussed in Ref. [13] the corrected neutral Higgs vertex can effect B s — B s mixing via 
double Higgs penguin diagrams that contribute to the Wilson coefficients of the operators 

0^ R ={b a P L s a ){J^P R s^), Of LL ={b a P L s a )(b p P L s p ), 0* RR =(b a P R s a )(bPp R s p ). (14) 

(for details of the operator basis we use see Ref. [31]). In the large tan/5 regime the dom- 
inant contribution arises from the Wilson coefficient C^, as the pseudoscalar and scalar 
contributions approximately cancel for both Cf LL and Cf RR . For non-zero 8ll and Slr the 
contributions to C% R are typically suppressed by a factor of the strange quark mass. For 
the insertions 5 RR and S R l, however, it is possible to avoid this suppression factor via the 
diagram where one Higgs penguin is mediated by chargino exchange and the other by gluino 
exchange. In the case of non-zero Srr, for example, C% R is given by 

/^ij r 8<y s AjY?Tni tan 4 /3 . . . , . . . . 

where k contains the effects of resuming tan f3 enhanced contributions 

K = \ fig) 

(l + e s tan/5) x [l + (e s + e y y; 2 )tan/?] 3 ' v ; 



3 Numerical Results 

When performing our numerical calculations we diagonalise all the various mixing matrices 
numerically, whilst we use FeynHiggs 2.0.2 [37] to compute the parameters associated with the 
Higgs sector. 4 Fig. 1 shows the effects of including beyond leading order GFM contributions 
on the decay B s — > and its dependence on the flavour violating parameters 5rr and 

5rl- For Srr the effect of using the bare bottom quark mass, as well as the additional effects 
induced by the off-diagonal elements of m$ present in Eq. (10), reduce the contribution 
to the decay by about a factor of two. A similar reduction between the LO and BLO 
results occurs for the insertion S LL . This reduction of LO effects by the inclusion of BLO 

4 We therefore ignore any possible effects on the Higgs sector due to off-diagonal entries in the squark 
mass matrix. However, the corrections due to 5ll to the lightest Higgs mass tend to be rather small [38]. 
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Figure 1: BR(i? s — > ) vs. (on the left) and vs. Srl (on the right). The soft sector is 
parameterised in terms of a common mass for the squark soft terms rriq, m~ g = \/2mq = 1 TeV, 



A t = -ITeV, m H + 
tan/5 = 50. 



jj, = 500 GeV and the gaugino soft terms Mi = M 2 = 0.5m,g, for 



supersymmetric corrections can be viewed as an extension of the focusing effect found in 
Ref. [14] to the neutral Higgs vertex. For 5rl (and similarly 5lr) the effects are more 
dramatic. As stated above, at leading order the contribution due to 5rl cancels and the 
remaining contributions stem from Z penguin diagrams that are not enhanced at large tan f3. 
However, once BLO corrections are taken into account the effects can differ quite significantly 
from the LO scenario where the dominant contributions at large tan /3 arise from the chargino 
contributions to the neutral Higgs vertex. 




5 RR 8 RL 



Figure 2: AM Bs vs. 5rr (on the left) and vs. Srl (on the right) for the same parameters as 
in Fig. 1. 
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Fig. 2 shows a similar plot for the B s — B s system. The strong linear dependence on both 
5rl and 8 RR confirms the approximate formula presented in the previous section. The effect 
of including BLO contributions once again has a large effect. Both graphs are somewhat 
similar to their B s — > analogs underlining the dependence on the corrected Higgs 

vertex and the large effects it can display in the large tan/3 region. For 8 RR the effect of 
including BLO contributions can once again lessen the dependence on 8 RR with respect to a 
purely LO analysis and lead to values of AMb s closer to the SM prediction. In the case of 
the insertion 8 RL , the effect of including BLO contributions is once again rather large and 
can lead to large deviations from a purely LO calculation. 

From Figs. 1-2 it is evident that the inclusion of BLO effects can significantly reduce 
both BK(B S — > and AMg s relative to LO predictions in the case of the insertion 

5 RR . This focusing effect can be viewed as an extension of the results presented in Ref. [14] 
to the neutral Higgs vertex. Fig. 3 illustrates the strong correlation between AM Bg and 
BR(B S — > at large tan (3 reflecting the fact that both processes are highly dependent 

on the neutral Higgs coupling in this regime. 
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Figure 3: A scatter plot of BR(B S -> vs. AM Bs for 0.2 < \5 RR \ < 0.5 and tan/3 = 40. 

Finally let us consider the effects the current experimental bounds on B s — > fi + fi~ and 
AM Bs have on the flavour violating parameters. As discussed in Ref. [14] the constraints 
placed on 8 RR and 8 R l by B — > X s 7 are rather weak. This is mainly due to the fact that the 
dominant contributions from both insertions are to the primed Wilson coefficients C 78 . Since 
there is no interference between the primed and unprimed operators the contributions to the 
final branching ratio are always positive. These additional contributions can therefore, in 
the CMSSM (mSUGRA) favoured scenario A t < 0, fi > 0, counter the chargino contribution 
which tends to decrease the overall branching ratio. These arguments however do not apply 
to the B s — B s system and the decay B s — > where it is quite possible for the behaviour 

at large tan (3 to be completely dominated by the effects of the Higgs penguin contributions. 
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Figure 4: In the above plots the region excluded by B — > is shaded in yellow (light 
grey). The additional regions excluded by B s — > /i + ^ and AMb 3 are shaded in orange 
(medium grey) and red (dark grey), respectively. The soft sector is parameterised as follows, 
rrig = V2rriq, A t = —nig and [i = m H + = 0.5m^, for tan/3 = 40. 

In Fig. 4 we show how the additional constraints supplied by the decay B s — > and 
AMB a affect the otherwise permitted values of both 5rl and 5rr. In both plots we used the 
current 95% confidence limits on AMb 3 and B s — > [i + , 

BR(S S -> < 5.0 x 10~ 7 , AM Bs > 14.5ps -1 . (17) 

For B — > X s 7 we combined the current experimental and theoretical errors in quadrature 
and added a small additional error to account for the accuracy of the supersymmetric portion 
of the calculation 

2.72 x 10~ 4 < BR(B -> X s7 ) < 3.96 x 10" 4 . (18) 

It can be seen from both plots that B s — > [i + fi~ and AM Bg can provide sizable constraints 
compared to B — > Xg'j. AMb s in particular provides a rather effective means of constraining 
5rl and Srr. As stated in the previous section, the contributions due to 8rl and 5rr to 
the Wilson coefficient C% R are not suppressed by factors of m s . Coupled with the tan 4 /3 
dependence of the contribution the effects of the double penguin can be rather large and can 
compete with the Standard Model contribution at large tan j3. 

4 Conclusions 

We have found that by taking into account GFM contributions when calculating the ra- 
diative corrections to the down quark mass matrix, the 5lr (and 5rl) dependence of the 
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corrected neutral Higgs vertex that conventionally cancels in LO calculations can reappear. 
The behaviour of processes that are highly dependent on this vertex (such as B s — > ji + /i~ 
and B s — B s mixing) can therefore change dramatically once GFM corrections are taken into 
account. In the case of the insertions 8rr and 5ll the effects of including beyond leading or- 
der GFM contributions typically reduce the values of BK(B S — > fi + fj>~) and AMb s compared 
to a purely LO calculation, exhibiting a focusing effect in the Higgs sector similar to the one 
pointed out in the case of B — > X s ^ in Ref. [14]. 

In the second part of our analysis we have illustrated how these effects can constrain 
the values of the flavour violating parameters 5rl and 5rr- The strong enhancement that 
super symmetric contributions to AM# s and BK(B S — » receive for non-zero 8rl and 

5rr, can lead to far stricter constraints on these parameters than in an analysis that just 
takes into account the effects that they have on B — > X s ~f. 
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A Loop Functions 

The loop functions H 2 (xi,x 2 ) and H 3 (xi,x 2 ,x 3 ) are given by 

H j Xi x \ x 1 \ogx 1 + x 2 hgx 2 , 1Q s 

2 U 2 (1 - Xi) (xi - x 2 ) (1 - x 2 ) (x 2 - £1) 
„ / x H 2 (x u x 2 ) - H 2 (x u x 3 ) 

H 3 (x 1: x 2 ,x 3 ) = (20) 

x 2 - x 3 
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